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Abstract 

A  simple  model  based  on  the  mean  spherical  approximation  is  described  for  the  dielectric 
properties  of  polar  solvents.  The  molecular  properties  required  are  the  dipole  moment,  the 
molecular  polarizability,  and  the  molecular  diameter.  It  is  shown  that  inclusion  of  a  stickiness 
parameter  accounting  for  other  electrostatic  features  and  non-sphericity  in  the  real  system  leads  to  a 
simple  analytical  result  from  which  the  dielectric  constant  may  be  calculated.  The  stickiness 
parameter  is  shown  to  be  correlated  to  the  magnitude  of  dipole-dipole  interactions  in  the  system.  It 
either  aligns  or  misaligns  the  electrical  dipoles  of  the  solvent,  and  therefore  can  be  called  a  stnicture 
forming  or  struciute  breaking  parameter. 
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Introduction 

Application  of  the  mean  spherical  approximation  (MSA)  to  the  description  of  the  dielecu  ic 
properties  of  polar  solvents  dates  from  the  work  of  Wertheim.  ‘  In  this  nxxlcl,  solvent  molecules 
are  represented  as  hard  spheres  with  centrally  located  point  dipoles.  When  applied  to  water,  the 
estimated  dielectric  constant  is  considerably  less  than  the  experimental  valuc.^  Thus,  one  must 
replace  the  experimental  dipole  moment  by  a  higher  effective  value  in  order  to  obtain  ilie 
experimental  dielectric  constant  from  the  MSA  model.23  It  is  generally  recognized  that  the  failure 
of  models  for  water  based  on  a  hard  sphere  -  pioint  dipiole  description  is  due  to  the  fact  tliat  this 
simple  molecule  has  a  more  complex  multipxjle  character.'*  The  complex  electrostatic  dcscnptjon 
leads  to  its  well  known  hydrogen-bonded  character. 

Application  of  the  MSA  model  to  larger  pxilar  solvents  is  certainly  expected  to  result  in  other 
problems.  Thus,  if  one  considers  the  series  of  normal  alcohols,  the  representation  of  the  molecule 
as  a  sphere  with  a  central  point  dipole  becomes  worse  as  the  length  of  the  alkyl  chain  incrca.scs 
Problems  of  a  similar  nature  arc  encountered  when  one  considers  a  scries  of  aptotic  solvents  .such 
as  acetonitrile,  and  alkyl  cyanides  with  larger  alkyl  groups.  Thus,  in  general,  one  recognizes  two 
problems  with  any  model  for  a  pxrlar  liquid  based  on  hard  spheres  with  centrally  located  p>oint 
dijxjles.  The  first  is  that  the  electrostatic  description  of  the  molecule  may  require  that  the 
quadrupole  moment  or  higher  order  multiplies  be  considered.  The  second  problem  is  that  the 
molecule  may  not  be  well  represented  as  a  sphere. 

One  way  of  dealing  with  the  above  inadequacies  is  to  introduce  stickiness  in  the  model.  In 
this  way,  one  can  correct  for  the  multipole  nature  of  the  pwlar  molecule  and  its  soft  iniermolecular 
interactions.  The  sticky  interaction  either  aligns  or  misaligns  the  electrical  dipx)les  of  the  solvent 
molecules,  and  therefore  can  be  called  a  structure  forming  or  structure  brc<tking  parameter.  In  the 
present  paper,  the  original  work  of  Wertheim*  is  extended  by  including  a  sticky  interaction 
parameter  in  the  description  of  the  polar  molecule.  Otherwise,  the  solvent  molecule  is  considered 
to  be  a  hard  sphere  possessing  a  permanent  point  dipole  and  a  polarizability.  The  model  is  applied 


to  representative  polar  solvents,  including  both  aproric  and  prone  systems,  in  order  to  assess  thc 
role  of  the  stickiness  parameter  in  determining  bulk  dielectric  properties. 

The  Model 

On  the  basis  of  Wenheim's  original  work*,  one  defines  a  polarization  parartKicr,  A,  by  the 
relationship 

X2(X+ 1)4=  I6es  (1) 

where  £§  is  the  static  dielectric  constant  of  the  polar  liquid.  The  parameter,  X,  also  arises  when  the 
MSA  is  applied  to  the  estimation  of  ion-dipole  interactions  in  an  infinitely  dilute  electrolyte 
solution^,  and  is  designed  to  take  into  account  softness  and  nonspherccity  of  die  solvent  core,  a  is 
obtained  from  the  excess  energy  parameter,  b2,  which  comes  from  the  solution  of  the  MSA  for 
dipole-dipole  interactions  for  a  system  composed  of  hard  spheres  with  centrally  located  point 
dipoles.  The  relationship  between  X  and  b2  is 


1  _  l~^b2/3 
^  “  l-b2/6 


(2) 


where 

PO 

b2  =  37tps-\/^  J  dr^r^ 

o 

h^^2  (j.)  is  the  relevant  coefficient  from  the  invariant  expansion  of  the  total  pair  correlation  function 
h{r)  which  depends  on  the  distance,  r,  between  two  dipolar  spheres.^  Ps  is  the  number  density  of 
molecules  in  the  liquid  and  Os,  their  diameter.  The  excess  internal  energy  per  unit  volume  for  the 
liquid  is  given  by 


E 

V  kBT  - 


d^b2 

2n  Os^ 


(4) 


where 


47t  Ps  p2 


Here,  p,  is  the  permanent  dipole  moment  of  the  solvent  nwlecule,  and  kjj,  the  Boltzmann  constant 
In  the  MSA,  when  the  molecules  are  non-polarizable,  one  obtains  the  result 


(X+2)2 
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(6) 


Wertheim  also  studied  the  MSA  and  other  approximations  for  the  case  of  dipole  hard  spheres 
which  are  polarizable.^  This  work  was  reformulated  to  show  that  it  is  equivalent  to  a  mixture  of 
hard  spheres  of  equal  diameter  but  varying  dipole  moment. 2  A.s  a  result,  when  the  molecu]e.s  have 
a  polarizability  a,  eq.  (6)  becomes 


where 
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with 


-1  = 


4  b2  g 
30s3 


(9) 


Solving  eq.  (7)  for  the  parameter  X,  one  obtains 

^  =  (1  +  3da)'^  -  1  (10) 

where 

da  =  dp/(l-l/es)»/2  (11) 

When  one  applies  the  above  system  of  equations  to  water  for  which  p  is  equal  to  1 .83 
Debyes  and  a,  to  1.47x10'^  nm^,  one  obtains  the  result  that  dp  =  4.33  and  X  =  2.75.  As  a  result, 
on  the  basis  of  eq.  (1),  the  MSA  estimate  of  Es  for  water  at  25  °C  is  93.7,  which  can  be  compared 
with  the  experimental  result  of  78.3.  This  suggests  that  a  small  correction  to  the  contact  correlation 
functions  would  yield  good  agreement  with  experiment.  This  is  achieved  here  using  the  directional 
sticky  model  described  by  Blum,  Cummings  and  Bratko.8  In  its  simplest  dipolar  version. 


introducing  a  longitudinal  stickiness  parameter,  to,  and  a  transversal  suckiness  parameter,  tj.  eq 
(7)  becomes^ 
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where  P3,  ^6  and  Pn  are  functions  of  X  namely 


P3  = 

P6  = 


3X 

X+2 

3 

X+2 


(13) 

(14) 


and 


P12  = 


3(X+1) 

2(U2) 


(15) 


Since  the  transversal  interactions  do  not  contribute  to  the  dielectric  constant,  ti,  may  be  set  equal  to 
zero.  Substituting  for  the  P’s  on  the  basis  of  eqs.  (13)-(15),  eq.  (12)  becomes 
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(16) 


or,  using  eq.  (1), 

07) 

By  adding  the  longitudinal  stickiness  parameter  to  the  estimation  of  dipole-dipole  interactions,  we 
provide  a  simple  and  easy  way  to  correct  for  the  deficiencies  of  the  MSA  in  the  close  range  region 
of  the  pair  correlation  function.  In  addition,  other  non-spherical  forces  effecting  dielectric 
behaviour  are  dealt  with. 

Having  estimated  dp  on  the  basis  of  eq.  (8),  one  can  then  determine  the  stickiness  parameter 
to  by  solving  eq.  (17)  using  the  value  of  X  obtained  from  the  experimental  value  of  £5  (eq.  (D). 
One  then  has  a  self-consistent  model  in  which  the  experimental  value  of  £5  is  estimated  from  the 


molecular  dipole  moment  and  polarizability  using  the  experimentally  defined  polarization  paranveicr 
X  with  an  additional  stickiness  parameter  t<j. 

This  model  is  now  examined  for  selected  polar  solvents.  However,  before  examining  the 
static  solvent  permittivity,  the  high  frequency  permittivity,  and  the  ability  of  the  MSA  to  estimate 
the  refractive  index  firom  the  molecular  polarizability  arc  considered. 

Solvent  Permittivity  at  Optical  Frequencies 

At  optical  frequencies,  the  system  only  responds  to  an  electrical  field  via  its  molecular 
polarizability.  Thus,  in  this  case  the  parameter  dp  is  given  by  the  equation 

dp  =  4jta  pstt  (IH) 

Furthermore,  the  polarization  parameter,  X,  is  now  calculated  from  the  optical  permittivity,  £op, 
and  thus  becomes  Xop: 

eoF  =  (19) 

where  Eop  is  the  square  of  the  refractive  index.  Thus,  an  alternate  way  to  estimating  dp  is  via  the 
equation 

^^(Xc,p+2)2 

a,  . 

Since  permanent  dipole-dipole  interactions  do  not  contribute  to  the  dielectric  properties  at  high 
frequencies,  the  stickiness  parameter  to  does  not  appear  in  eq.  (20).  The  parameters  bx  and  a 
defined  above  by  equations  (2)  and  (9),  respectively,  arc  now  estimated  using  the  value  of  Xop 
obtained  from  eq.  (19). 

In  the  case  of  water  at  25  ®C,  the  value  of  eop  is  1 .7756.  The  corresponding  value  of  the 
polarization  parameter  Xop  is  1.151  on  the  basis  of  eq.  (19).  One  may  now  estimate  the  parameter 
dp  by  using  eq.  (20)  and  then  calculate  the  molecular  polarizability  on  the  basis  of  eq.  (18).  For 
water,  the  MSA  estimate  of  a  is  1.483  x  lO'^  nm^,  a  result  which  compares  very  well  with  the 
value,  1.471  x  10'^  nm^,  estimated  on  the  basis  of  the  Lorentz-Lorenz  relationship.  ^  *  In  the  latter 
case,  the  polarizability  is  calculated  firom  the  molar  refraction  R  defined  by  the  equation 


R  =  ^  9  =  :|  n  No  a 
eop+2  ^ 


(21) 


where  V  is  the  molar  volume  and  Nq,  the  Avogadro  constant.  Keeping  in  mind  the  fact  that  tfic 
value  of  the  polari2ability  depends  very  much  on  whether  it  is  obtained  from  gas  or  liquid  phase 
measurements it  seems  appropriate  that  it  be  obtained  from  the  refractive  index  of  the  liquid. 

Values  of  a  were  calculated  for  all  other  solvents  considered  using  both  the  MSA  and  the 
Lorentz-Lx)renz  relationship  from  tabulated  values  of  Eop  and  In  all  cases,  the  MSA  esiiniate 

is  slightly  higher,  typically,  by  approximately  \%.  Because  the  difference  is  so  small,  the  MSA 
result  was  used  in  all  further  calculations  and  is  tabulated  in  Table  1  together  with  other  pertinent 
solvent  data. 

The  Static  Solvent  Permittivity 

The  present  model  was  applied  to  20  different  polar  solvents,  both  proiic  and  aproiic,  chosen 
so  that  the  role  of  the  stickiness  parameter  could  be  assessed.  Three  solvent  groups  considered 
were  water  and  the  lower  alcohols  in  which  the  pwlar  group  is  -OH,  the  alkyl  cyanides  with  -C  =  N 
as  the  polar  group,  and  the  amides  whose  dielectric  parameters  depend  markedly  on  whether  the 
solvent  is  protic  or  aprotic.  The  dielectric  and  other  relevant  parameters  for  these  systems  are 
summarized  in  Table  1.  Some  comment  should  be  made  about  the  values  of  the  dipole  moment 
chosen  from  the  compilations  of  McClellan. The  value  of  this  parameter  depends  somewhat  on 
whether  it  was  measured  in  the  gas  phase  or  in  a  dilute  solution  in  a  non-polar  solvent.  The  chosen 
values  for  a  group  of  solvents  were  obtained  by  one  method,  preferably  using  measurements  from 
one  laboratory.  The  results  obtained  in  the  present  study  depend  markedly  on  the  value  chosen  for 
this  parameter,  so  it  is  important  that  it  properly  reflect  trends  within  a  given  group  of  solvents. 

The  stickiness  parameter,  to,  also  recorded  in  Table  1  was  estimated  from  the  dielectric  data 
in  the  following  way.  Initially,  the  polarization  parameter,  X,  and  the  related  energy  parameter,  b2, 
were  calculated  from  the  bulk  dielectric  constant,  Cs,  (eqs.  (1)  and  (2)).  Then  the  molecular 
parameter,  d2,  was  estimated  on  the  basis  of  eqs.  (5),  (8)  and  (9)  using  the  molecular  polarizability 


X 


found  earlier,  and  the  dipole  moment,  density,  and  molecular  diameter  tabulated  in  Table  1 
Finally,  the  value  of  to  was  obtained  by  solving  eq.  (17)  using  the  above  values  of  X  and  dp. 

From  the  results,  it  is  apparent  that  to  varies  over  a  wide  range  from  a  low  of  - 1 .498  for 
formamide  to  a  high  of  0.495  for  1-propanol.  Negative  values  of  to  indicate  that  dp  is 
underestimated  when  it  is  calculated  on  the  basis  of  eq.  (8).  It  also  follows  that  the  MSA  then 
underestimates  the  ordering  effect  of  a  given  dipole  on  its  neighbours.  In  the  case  of  proiic 
solvents  such  as  water  and  formamide,  this  is  attributed  to  the  presence  of  hydrogen  bonding,  and 
important  multipole  contributions  to  the  electrostatic  description  of  the  molecule.  Negative  values 
of  to  are  also  found  for  all  aptotic  solvents  considered.  This  also  reflects  the  fact  that  the  dipole 
moment  cannot  be  considered  to  be  centrally  located  in  any  of  the  polar  solvent  molecules 
considered. 

It  is  interesting  to  follow  the  change  in  to  with  the  strength  of  the  dipole-dipole  interactions  in 
the  system.  This  can  be  seen  for  the  water-alcohol  series  and  also  the  alkyl  nitrile  senes.  Plots  of 
to  against  p^/d^  are  linear  for  both  series  over  a  range  of  4  kT  (Fig.  1).  What  is  interesting  about 
these  results  is  that  to  becomes  more  positive  as  dipole-dipole  interactions  decrease  in  magnitude. 
This,  in  turn,  can  be  related  to  the  dipolar  group  in  the  solvent  molecule  moving  further  and  funher 
away  from  the  molecular  center.  No  such  trend  can  be  seen  for  the  amides.  Undoubtedly,  this  is 
due  to  the  fact  that  the  group  of  solvents  considered  includes  botl  protic  and  aptotic  systems  so 
that  the  complete  electrostatic  description  varies  considerably  among  the  members  of  this  group. 
N-methylformamide  is  the  only  other  solvent  listed  which  has  a  positive  value  of  t©.  Both 
formamide  and  N-methylformamide  are  systems  which  are  extensively  hydrogen  bonded  but  wiili 
quite  different  structures  as  a  result  of  change  in  number  of  acidic  hydrogens.  This  change 
undoubtedly  leads  to  the  large  difference  in  to  for  these  solvents. 

The  remaining  solvents  listed  in  Table  1  are  all  aptotic  and  possess  negative  values  of  to- 
This  parameter  is  plotted  against  the  dipole-dipole  interaction  energy  unit,  p^/d^  in  Fig.  2.  There  is 
clearly  a  correlation  between  these  two  quantities,  the  value  of  to  becoming  more  negative  as  pVd^ 
increases.  This  suggests  that  the  smallest  molecules  with  largest  dipole  moments  depart  most 


severely  from  the  simple  MSA  result.  This  result  reflects  the  fact  that  the  tendency  for  k>cal  order 
is  underestimated  by  the  MSA.  The  scatter  is  much  larger  tlian  that  seen  in  Fig.  1  where  molecules 
with  the  same  polar  group  were  considered.  Not  only  are  the  polar  groups  different  but  the  size 
and  shape  of  the  non-polar  part  of  the  molecule  varies  considerably  among  the  solvents  considered. 
In  addition,  systems  like  nitrobenzene  and  benzonitrile  have  very  polarizable  phenyl  rings  wherciis 
the  other  molecules  do  not.  In  conclusion,  the  observed  correlation  is  chiefly  a  result  of  the  change 
in  molecular  size,  and  corrects  the  MSA  model  for  the  fact  that  the  dipolar  group  is  not  centrally 
located  in  the  solvent  molecule. 

Discussion 

The  MSA  model  for  the  dielectric  properties  of  polar  solvents  described  here  is  remarkably 
simple.  Using  the  molecular  dipole  moment,  polarizability,  and  diameter  estimated  from  gas 
solubilities,  one  is  able  to  calculate  the  relative  permittivity  with  two  additional  parameters,  namely, 
the  polarization  parameter,  X,  and  the  stickiness  parameter,  to-  The  first  of  these,  X,  is  calculated 
direedy  from  the  static  dielectric  constant.  Furthermore,  the  stickiness  parameter,  to,  is  clearly 
correlated  to  the  ratio  pVd^  which  is  proportional  to  the  dipole-dipole  interaction  energy.  As  the 
size  of  the  non-polar  group  increases  in  a  series  of  molecules  with  the  same  p>olar  group,  the  ratio 
p2/d3  decreases,  and  to,  increases.  The  correlation  between  these  quantities  is  perfectly  linear  for  a 
given  series  of  compounds. 

The  limitations  of  the  simple  MSA  in  describing  the  properties  of  polar  liquids  are  well 
known  the  present  extension  adding  very  little  to  the  complexity  of  the  model.  If  one  compares 
the  results  of  the  present  model  to  those  from  molecular  dynamics  calculations,  it  is  clear  that  the 
latter  approach  only  would  yield  reasonable  results  if  a  detailed  electrostatic  description  of  the 
solvent  molecule  is  provided.  The  fact  that  the  MSA  provides  an  analytical  model  with  empirical 
parameters  which  can  be  qualitatively  related  to  known  properties  of  the  system  suggests  that  it  is  a 
very  good  basis  for  describing  the  dielectric  properties  of  polar  liquids. 
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Table  1.  Summary  of  Relevant  Solvent  Parameters 
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Legends  for  Figures 

Figure  1.  Plot  of  the  stickiness  parameter,  to,  against  the  dipole-dipole  imcraciion  energy  factor 
p^/d^  for  the  lower  alcohols  and  water,  and  for  the  alkyl  cyanides.  The  abbreviations 
for  the  solvent  given  next  to  each  point  are  defined  in  Table  1. 

Figure  2.  As  in  Figure  I,  but  for  all  aprotic  solvents  listed  in  Table  1. 


